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Abstract 

In this paper we propose the time-dependent generalization of an 'ordinary' autonomous 
human biomechanics, in which total mechanical + biochemical energy is not conserved. We 
introduce a general framework for time-dependent biomechanics in terms of jet manifolds 
associated to the extended musculo-skeletal configuration manifold, called the configuration 
bundle. We start with an ordinary configuration manifold of human body motion, given as a 
set of its all active degrees of freedom (DOF) for a particular movement. This is a Riemannian 
manifold with a material metric tensor given by the total mass-inertia matrix of the human 
body segments. This is the base manifold for standard autonomous biomechanics. To make its 
time-dependent generalization, we need to extend it with a real time axis. By this extension, 
using techniques from fibre bundles, we defined the biomechanical configuration bundle. On 
the biomechanical bundle we define vector-fields, differential forms and afline connections, as 
well as the associated jet manifolds. Using the formalism of jet manifolds of velocities and 
accelerations, we develop the time-dependent Lagrangian biomechanics. Its underlying geo- 
metric evolution is given by the Ricci flow equation. 

Keywords: Human time-dependent biomechanics, configuration bundle, jet spaces, Ricci flow 

1 Introduction 

It is a well-known fact that most of dynamics in both classical and quantum physics is based 
on assumption of a total energy conservation (see, e.g. [I]). Dynamics based on this assump- 
tion of time-independent energy, usually given by Lagrangian or Hamiltonian energy function, is 
called autonomous. This basic assumption is naturally inherited in human biomechanics, formally 
developed using Newton-Euler, Lagrangian or Hamiltonian formalisms (see El HI El El 13 [HI E] ) - 
And this works fine for most individual movement simulations and predictions, in which the 
total human energy dissipations are insignificant. However, if we analyze a lOOm-dash sprinting 
motion, which is in case of top athletes finished under 10 s, we can recognize a significant slow-down 
after about 70 m in all athletes despite of their strong intention to finish and win the race, which 
is an obvious sign of the total energy dissipation. This can be seen, for example, in a current 
record-braking speeddistance curve of Usain Bolt, the world-record holder with 9.69 s [10], or in 
a former record-braking speeddistance curve of Carl Lewis, the former world-record holder (and 9 
time Olympic gold medalist) with 9.86 s (see Figure 3.7 in [9]). In other words, the total mechanical 



+ biochemical energy of a sprinter cannot be conserved even for 10 s. So, if we want to develop a 
realistic model of intensive human motion that is longer than 7-8 s (not to speak for instance of a 
4 hour tennis match), we necessarily need to use the more advanced formalism of time-dependent 
mechanics. 

Similarly, if we analyze individual movements of gymnasts, we can clearly see that the high 
speed of these movements is based on quickly- varying mass-inertia distribution of various body 
segments (mostly arms and legs). Similar is the case of pirouettes in ice skating. As the total mass- 
inertia matrix My of a biomechanical system corresponds to the Riemannian metric tensor g^ of 
its configuration manifold, we can formulate this problem in terms of time-dependent Riemannian 
geometry [Ud]. 

The purpose of this paper is to introduce a general framework for time-dependent biomechanics, 
consisting of the following four steps: 

1. Human biomechanical configuration manifold and its (co)tangent bundles; 

2. Biomechanical configuration bundle, as a time-extension of the configuration manifold; 

3. Biomechanical jet spaces and prolongation of locomotion vector-fields developed on the con- 
figuration bundle; and 

4. Time-dependent Lagrangian dynamics using biomechanical jet spaces. 

In addition, we will show that Riemannian geometrical basis of this framework is defined by the 
Ricci flow. In particular, we will show that the exponential-like decay of total biomechanical 
energy (due to exhaustion of biochemical resources [9] ) is closely related to the Ricci flow on the 
configuration manifold of human motion. 

2 Configuration Manifold for Autonomous Biomechanics 

Recall from [5] that representation of an ideal humanoid-robot motion is rigorously defined in 
terms of rotational constrained S'0(3)-groups in all main robot joints. Therefore, the configuration 
manifold Q ro b for humanoid dynamics is defined as a topological product of all included 50(3) 
groups, Qrob = Yii SO(3y. Consequently, the natural stage for autonomous Lagrangian dynamics 
of robot motion is the tangent bundle TQ ro t,, defined as follows. To each n— dimensional (nD) 
configuration manifold Q there is associated its 2nD velocity phase-space manifold, denoted by 
TQ and called the tangent bundle of Q. The original smooth manifold Q is called the base of TQ. 
There is an onto map 7r : TQ — + Q, called the projection. Above each point x £ Q there is a 
tangent space T X Q = tt^ 1 (x) to Q at x, which is called a fibre. The fibre T X Q C TQ is the subset 
of TQ, such that the total tangent bundle, TQ = | | T X Q, is a disjoint union of tangent spaces 

m£<2 

T X Q to Q for all points x € Q. From dynamical perspective, the most important quantity in the 
tangent bundle concept is the smooth map v : Q — * TQ, which is an inverse to the projection tt, 
i.e, ttov = Idg, 7r(i>(x)) = x. It is called the velocity vector-field. Its graph (x,v(x)) represents the 
cross-section of the tangent bundle TQ. This explains the dynamical term velocity phase-space, 
given to the tangent bundle TQ of the manifold Q. The tangent bundle is where tangent vectors 
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live, and is itself a smooth manifold. Vector-fields are cross-sections of the tangent bundle. Robot's 
Lagrangian (energy function) is a natural energy function on the tangent bundle TQQ 

On the other hand, human joints are more flexible than robot joints. Namely, every rotation 
in all synovial human joints is followed by the corresponding micro-translation, which occurs after 
the rotational amplitude is reached [B]. So, representation of human motion is rigorously defined in 
terms of Euclidean S'i?(3)-groups of full rigid-body motion [51 SI [T] in all main human joints (see 
Figure [l}. Therefore, the configuration manifold Q for human dynamics is defined as a topological 
product of all included constrained S , £ , (3)-groups0 Q = ^SE^S) 1 . Consequently, the natural 
stage for autonomous Lagrangian dynamics of human motion is the tangent bundle TQ (and for 
the corresponding autonomous Hamiltonian dynamics is the cotangent bundle T*Q). 

Therefore, the configuration manifold Q for human musculo-skeletal dynamics is defined as a 
Cartesian product of all included constrained SE(3) groups, Q — JTj SE(3y where j labels the 
active joints. The configuration manifold Q is coordinated by local joint coordinates x % (t), i = 
f, ...,n = total number of active DOF. The corresponding joint velocities x l (t) live in the velocity 
phase space TQ, which is the tangent bundle of the configuration manifold Q. 

The velocity phase space TQ has the Riemannian geometry with the local metric form: 

(g) = ds 2 = gijdx' dx 1 , (Einstein's summation convention is in always use) 

where g%j(x) is the material metric tensor defined by the biomechanical system's mass-inertia 
matrix and dx l are differentials of the local joint coordinates x 1 on Q. Besides giving the local 
distances between the points on the manifold Q, the Riemannian metric form (g) defines the 
system's kinetic energy: 

T = —giji x 1 , 

1 The corresponding autonomous Hamiltonian robot dynamics takes place in the cotangent bundle T*Q ro i,, 
defined as follows. A dual notion to the tangent space T m Q to a smooth manifold Q at a point m is its cotangent space 
T^Q at the same point m. Similarly to the tangent bundle, for a smooth manifold Q of dimension n, its cotangent 
bundle T*Q is the disjoint union of all its cotangent spaces T^Q at all points in 6 Q, i.e., T*Q = | | T^Q. 

mGQ 

Therefore, the cotangent bundle of an n— manifold Q is the vector bundle T*Q = (TQ)*, the (real) dual of the 
tangent bundle TQ. The cotangent bundle is where 1-forms live, and is itself a smooth manifold. Covector— fields 
(1-forms) are cross-sections of the cotangent bundle. Robot's Hamiltonian is a natural energy function on the 
cotangent bundle. 

2 Briefly, the Euclidean SE(3)-group is defined as a semidirect (noncommutative) product (denoted by >) of 3D 
rotations and 3D translations: SE(3) := SO(3) >R 3 . Its most important subgroups are the following (for technical 
details see HimiTl): 



Subgroup 




Definition 




50(3), group of rotations 
in 3D (a spherical joint) 


Set of all proper orthogonal 
3 X 3 — rotational matrices 


SE(2), special Euclidean group 
in 2D (all planar motions) 


S< 


it of all 3 X 3 — matrict 
cos 9 sin 9 r x 
— sin 9 cos 9 r y 
1 


:s: 


SO(2), group of rotations in 2D 
subgroup of S£(2)-group 
(a revolute joint) 


Set of all proper orthogonal 
2 X 2 — rotational matrices 
included in SE(2) — group 


R 3 , group of translations in 3D 
(all spatial displacements) 


Euclidean 3D vector space 
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Figure 1: The configuration manifold Q of the human musculoskeletal dynamics is defined as 
a topological product of constrained SE(3) groups acting in all major (synovial) human joints, 
Q = UiSE(3)\ 

giving the Lagrangian equations of the conservative skeleton motion with kinetic-minus-potential 
energy Lagrangian L = T — V , with the corresponding geodesic form [5] 

j f L ± ,-L x ,=0 or r ■ l',.i-r h 0, (1) 

where subscripts denote partial derivatives, while Tj k are the Christoffel symbols of the affine 
Levi-Civita connection of the biomechanical manifold Q. 

This is the basic geometrical structure for autonomous Lagrangian biomechanics. In the next 
sections will extend this basic structure to embrace the time-dependent Lagrangian biomechanics. 

3 Biomechanical Bundle, Sections and Connections 

While standard autonomous Lagrangian biomechanics is developed on the configuration manifold 
X, the time-dependent biomechanics necessarily includes also the real time axis R, so we have 
an extended configuration manifold Rxl. Slightly more generally, the fundamental geometrical 
structure is the so-called configuration bundle n : X — > R. Time-dependent biomechanics is thus 
formally developed either on the extended configuration manifold K x X , or on the configuration 
bundle 7r : X — > R, using the concept of jets, which are based on the idea of higher-order tan- 
gency, or higher-order contact, at some designated point (i.e., certain joint) on a biomechanical 
configuration manifold X. 

In general, tangent and cotangent bundles, TM and T*M, of a smooth manifold M, are special 
cases of a more general geometrical object called fibre bundle, denoted tt : Y — > X, where the word 
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fiber V of a map tt : Y — > X is the preimage 7r _1 (a;) of an element x E -X". It is a space which 
locally looks like a product of two spaces (similarly as a manifold locally looks like Euclidean space), 
but may possess a different global structure. To get a visual intuition behind this fundamental 
geometrical concept, we can say that a fibre bundle Y is a homeomorphic generalization of a 
product space X x V (see Figure ^ , where X and V are called the base and the /i&re, respectively. 
tt : Y — > X is called the projection, Y x — tt^ 1 (x) denotes a fibre over a point a; of the base X, while 
the map / = tt^ 1 : X — > K defines the cross-section, producing the graph (x, /(as)) in the bundle 
y (e.g., in case of a tangent bundle, f — x represents a velocity vector-field). 




Figure 2: A sketch of a locally trivial fibre bundle Y X x V a,s & generalization of a product 
space X x V; left - main components; right - a few details (see text for explanation). 

More generally, a biomechanical configuration bundle, tt : Y — > X, is a locally trivial fibred (or, 
projection) manifold over the base X. It is endowed with an atlas of fibred bundle coordinates 
(x x ,y l ), where (x x ) are coordinates of X. 

All dynamical objects in time-dependent biomechanics (including vectors, tensors, differential 
forms and gauge potentials) are cross-sections of biomechanical bundles, representing generaliza- 
tions of graphs of continuous functions. 

An exterior differential form a of order p (or, a p—form a) on a base manifold A is a section 
p 

of the bundle A T*X — > X [19 . It has the following expression in local coordinates on X 

a = a\ 1 ...x p dx Xl A • • ■ A dx Xp (such that |a|=p), 

where summation is performed over all ordered collections (Ai, X p ). f2 p (A) is the vector space 
of p— forms on a biomechanical manifold X. In particular, the 1-forms are called the Pfaffian 
forms. 

The contraction J of any vector-field u = u^d^ and a p—form a — ax 1 ...x p dx 1 A • • • A dx Xp on 
a biomechanical manifold X is given in local coordinates on X by 

u\a = u tJl a tl \ 1 ...\ p _ 1 dx Xl A • • • A dx Xp -\ 

It satisfies the following relation 

u \(a A p) = u\a A (3 + (-l) |a| a A u\/3. 

The Lie derivative C u a of p—form a along a vector-field u is defined by Cartan's 'magic' formula 
(see HH]): 

C u a — u\da + d(u\a). 
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It satisfies the Leibnitz relation 

C u (a A 0) = C u a A /3 + a A C u /3. 
A linear connection f on a biomechanical bundle Y — » X is given in local coordinates on F by 

T = dx x ®[d x -T i jX (x)y j d i \. (2) 
An affine connection T on a biomechanical bundle Y — > X is given in local coordinates on F 



by 

r = dx x ® p A + (-rvory + r A (a))ai]. 

Clearly, a linear connection T is a special case of an affine connection T. 



4 Biomechanical Jets 

A pair of smooth manifold maps, /i,/2 : M N (see Figure are said to be k— tangent (or 
tangent of order k, or have a fcth order contact) at a point a; on a domain manifold M, denoted 
by /i ~ f 2 , iff 

/l(ar) = fzix) called — tangent, 
9 x fi(x) = d x f 2 [x), called 1 - tangent, 
d xx fi(x) = d xx f 2 {x), called 2 - tangent, 
etc. to the order k 

In this way defined k—tangency is an equivalence relation. 

fv 




Figure 3: An intuitive geometrical picture behind the k— jet concept, based on the idea of a higher- 
order tangency (or, higher-order contact). 

A k~jet (or, a jet of order fc), denoted by j£f, of a smooth map / : Q — > N at a point ieQ 
(see Figure [3]), is defined as an equivalence class of fc— tangent maps at a;, 



j x f 



N = {/' : /' is k — tangent to / at x}. 
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For example, consider a simple function f : X — > Y, x ^ y = f(x), mapping the X— axis into 
the Y — axis in R 2 . At a chosen point x G X we have: 
a 0— jet is a graph: (x, f(x)); 
a 1— jet is a triple: (x, f(x), f'(x)); 
a 2— jet is a quadruple: (a;, f(x), f'(x), f"(x)), 
and so on, up to the order k (where f'(x) = , etc). 

The set of all k— jets from j£f : X — » Y is called the fc— jet manifold J k (X, Y). 

Formally, given a biomechanical bundle Y — > X, its first order jet manifold J l Y comprises the 
set of equivalence classes j^s, x e X, of sections s : X — > Y so that sections s and s' belong to the 
same class iff 

Ts \ Tm x= Ts' \ Tx x ■ 

Intuitively, sections s,s' 6 j'^s are identified by their values s*(x) = s' l (x) and the values of their 
partial derivatives d )Jl s l (x) = d )1 s' t (x) at the point x of X. There are the natural fibrations [TH] 

7Ti : J X Y 3 ils k-> x S X, 7T i : J^Y 3 i-> s(a;) <G Y. 

Given bundle coordinates (x x ,y l ) of Y, the associated jet manifold J X Y is endowed with the 
adapted coordinates 

ox' 

In particular, given the biomechanical configuration bundle Q — > R over the time axis R, the 
1—jet space J 1 (R, Q) is the set of equivalence classes j\s of sections s l : R — * Q of the configuration 
bundle Q — > R, which arc identified by their values s l (t), as well as by the values of their partial 
derivatives <9 f s l = d t s l (t) at time points ( £ 1. The 1—jet manifold J 1 (R, Q) is coordinated by 
{t,q l ,q l ), that is by (time, coordinates and velocities) at every active human joint, so the 1-jets 
are local joint coordinate maps 

fia-.R^Q, t^(t,q\q l ). 

The repeated jet manifold J 1 J 1 Y is defined to be the jet manifold of the bundle J 1 Y — > X. It 
is endowed with the adapted coordinates (x A , y\ y\, y^yVx^)- 

The second order jet manifold J 2 Y of a bundle Y — > X is the subbundle of J 2 Y — > J 1 !" defined 
by the coordinate conditions y A = j/* iA . It has the local coordinates (x x , y l , y\, y\ <fl ) together 
with the transition functions [19] 

dx a 

The second order jet manifold J 2 Y of Y comprises the equivalence classes j 2 s of sections s of 
Y —y X such that 

y\(f x s) = d x s\x), yl^s) = d^(x). 

In other words, two sections s, s' € j 2 s are identified by their values and the values of their first 
and second order derivatives at the point x G X. 
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In particular, given the biomechanical configuration bundle Q — > R over the time axis R, the 
2— je£ space J 2 (R, Q) is the set of equivalence classes jfs of sections s l : R — > Q of the configuration 
bundle 7r : Q — > R, which are identified by their values s z (t), as well as the values of their first and 
second partial derivatives, d t s l = d t s l (t) and d u s l = d t ts % (t), respectively, at time points t G R. 
The 2-jet manifold J 2 (R, Q) is coordinated by (t, q l , q l ,q l ), that is by (time, coordinates, velocities 
and accelerations) at every active human joint, so the 2— jets are local joint coordinate map^| 

jfaiR^Q, t^(t,q\q\q l ). 

5 Lagrangian Time— Dependent Biomechanics 
5.1 Jet Dynamics and Quadratic Equations 

The general form of time-dependent Lagrangian biomechanics with time- dependent Lagrangian 
function L(t;q l ;q l ) defined on the jet space X — J 1 (R, Q) = R x TQ, with local canonical coor- 
dinates: (t;q l ;q l ) = (time, coordinates and velocities) in active local joints, can be formulated as 

mm 

—Lqi-L q i=J : i (t,q ) q), (i = l,...,n), (3) 

where the coordinate and velocity partial derivatives of the Lagrangian are respectively denoted 
by L q i and Lqi. 

The most interesting instances of (|3]) are quadratic biomechanical equations, of the general 
form 

C=q 4 = a) k (qn<i 3 q k + + f (<f )• (4) 

They are coordinate-independent due to the affine transformation law of coordinates q 1 . Then, it 
is clear that the corresponding dynamical connection L^ is affine [H Q] : 

I> = dq a ® [d a + (F A0 («f ) + Ti^q^df]. 

This connection is symmetric iff T l Xt2 = . 

There is 1-1 correspondence between the affine connections L on the affine jet bundle 
J 1 (R, Q) — > Q and the linear connections K on the tangent bundle TQ — ► Q of the autonomous 
biomechanical manifold Q. This correspondence is given by the relation 

Any quadratic biomechanical equation (j4|) is equivalent to the geodesic equation [19] 

t = l, t = 0, g < = o5 fc (^)« i ^ + 6j.(^t + / < (0«, 
for the symmetric linear connection 

K = dq a <g>(d a + KUt,q l Wd^ 

3 For more technical details on jet spaces with their physical applications, see 17 18 19, 20 ). 
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on the tangent bundle TQ — > Q, given by the components 

K av = 0, K 1 q = f\ K l j = A'/o = K j\ = a)k- 

Conversely, any linear connection K on the tangent bundle TQ — > Q defines the quadratic dynam- 
ical equation 

if = K 3 l k qiq k + (K l 3 + K/o)? + Kq'q, 
written with respect to a given reference frame (t, q 1 ) = q^ (see [TH] for technical details). 

5.2 Local Muscle— Joint Mechanics 

The right-hand side terms J~i(t, q, q) of ([3|) denote any type of external torques and forces, including 
excitation and contraction dynamics of muscular-actuators and rotational dynamics of hybrid 
robot actuators, as well as (nonlinear) dissipative joint torques and forces and external stochastic 
perturbation torques and forces. In particular, we have [2 [3]): 

1. Synovial joint mechanics, giving the first stabilizing effect to the conservative skeleton 
dynamics, is described by the (q, g)-form of the Rayleigh-Van der Pol's dissipation function 

1 n 

R = sE^) 2 ^ + 

i=l 

where on and (i i denote dissipation parameters. Its partial derivatives give rise to the viscous- 
damping torques and forces in the joints 

point = dR j dq i^ 

which are linear in q % and quadratic in q l . 

2. Muscular mechanics, giving the driving torques and forces j:™ uscle = j:m uscle (t, g, q) with 
(i = I, ... ,71) for human biomechanics, describes the internal excitation and contraction dynamics 
of equivalent muscular actuators [12) . 

(a) Excitation dynamics can be described by an impulse force-time relation 

pimp = F o^ _ e -t/T-*) if stimulation >0 
pimp = F o e ~t/r z if stimulation =0, 

where F® denote the maximal isometric muscular torques and forces, while r, denote the associated 
time characteristics of particular muscular actuators. This relation represents a solution of the 
Wilkie's muscular active-state element equation [T3] 

fi + Tn = TSA, fi(0) = 0, 0<5<1, 

where fj, — fi(t) represents the active state of the muscle, T denotes the element gain, A corresponds 
to the maximum tension the element can develop, and S = S(r) is the 'desired' active state as a 
function of the motor unit stimulus rate r. This is the basis for biomechanical force controller. 

(b) Contraction dynamics has classically been described by the Hill's hyperbolic force-velocity 
relation [T3] 

pHill = (jfh ~ ) 
(Sijtf + h) 



9 



where a, and bt denote the Hill's parameters, corresponding to the energy dissipated during the 
contraction and the phosphagenic energy conversion rate, respectively, while 5ij is the Kronecker's 
6— tensor. 

In this way, human biomechanics describes the excitation/contraction dynamics for the ith 
equivalent muscle-joint actuator, using the simple impulse-hyperbolic product relation 

F™ uscle (t,q,q) = Fl mp X Ff iU . 

Now, for the purpose of biomedical engineering and rehabilitation, human biomechanics has 
developed the so-called hybrid rotational actuator. It includes, along with muscular and viscous 
forces, the D.C. motor drives, as used in robotics p~5l [2] 

H 0h ° = *k(*) - Mk{t) - B k q k (t), with 
hik(t) + Rkik{t) + C k q k (t) = u k (t), 

where k — 1, . . . , n, i k (t) and u k (t) denote currents and voltages in the rotors of the drives, R k , l k 
and C k are resistances, inductances and capacitances in the rotors, respectively, while J k and B k 
correspond to inertia moments and viscous dampings of the drives, respectively. 

Finally, to make the model more realistic, we need to add some stochastic torques and forces 

m 

Tf° ch = B ij [q i {t),t] dW*(t), 

where Bij [q(t) , t] represents continuous stochastic diffusion fluctuations, and (t) is an N— variable 
Wiener process (i.e., generalized Brownian motion), with 

dW j (t) = W\t + dt) - W j (t), (for j = 1, . . . , TV). 



6 Time— Dependent Riemannian Geometry of Biomechanics 

As illustrated in the introduction, the mass-inertia matrix of human body, defining the Rieman- 
nian metric tensor gij(q) need not be time-constant, as in case of fast gymnastic movements and 
pirouettes in ice skating, which are based on quick variations of inertia moments and products 
constituting the material metric tensor gij(q). In particular, in the geodesic framework JI]), the 
(in)stability of the biomechanical joint and center-of-mass trajectories is the (in)stability of the 
geodesies, and it is completely determined by the curvature properties of the underlying manifold 
according to the Jacobi equation of geodesic deviation [H [T] 

D 2 J' pi W k dq™ 
as as ds 

whose solution J, usually called Jacobi variation field, locally measures the distance between nearby 
geodesies; D /ds stands for the covariant derivative along a geodesic and R l j km are the components 
of the Riemann curvature tensor. 

In general, the biomechanical metric tensor g^ is both time and joint dependent, gij — gijit, q). 
This time-dependent Riemannian geometry can be formalized in terms of the Ricci flow |21j , the 
nonlinear heat like evolution metric equation: 

dtgij = -Rij, (5) 
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for a time-dependent Ricmannian metric g = gij{t) on a smooth n— manifold Q with the Ricci 
curvature tensor Rij. This equation roughly says that we can deform any metric on the config- 
uration manifold Q by the negative of its curvature; after normalization, the final state of such 
deformation will be a metric with constant curvature. The negative sign in (O insures a kind of 
global volume exponential decay^ since the Ricci flow equation (O is a kind of nonlinear geometric 
generalization of the standard linear heat equation 

dtu = Au. 

In a suitable local coordinate system, the Ricci flow equation §5§ on a biomechanical configu- 
ration manifold Q has a nonlinear heat-type form, as follows. At any time t, we can choose local 
harmonic coordinates so that the coordinate functions are locally defined harmonic functions in 
the metric g(t). Then the Ricci flow takes the general form [T] 

d t gij = A Q g tj +G i3 (g,dg), where (6) 

is the Laplace-Beltrami operator of the configuration manifold Q and Gij(g,dg) is a lower-order 
term quadratic in g and its first order partial derivatives dg. From the analysis of nonlinear heat 
PDEs, one obtains existence and uniqueness of forward-time solutions to the Ricci flow on some 
time interval, starting at any smooth initial metric go on Q. 

The exponentially-decaying geometrical diffusion ([5]) is a formal description for pirouettes in 
ice skating and fast rotational movements in gymnastics. 

7 Conclusion 

We have presented the time-dependent generalization of an 'ordinary' autonomous human biome- 
chanics, in which total mechanical + biochemical energy is not conserved. We have introduced 
a general framework for time-dependent biomechanics in terms of jet manifolds associated to the 
extended musculo-skeletal configuration manifold, called the configuration bundle. We start with 
an ordinary, autonomous configuration manifold of human body motion, given as a set of its all 
active DOF for a particular movement. This is a Riemannian manifold with a material metric ten- 
sor given by the total mass-inertia matrix of the human body segments. This is the base manifold 
for standard autonomous biomechanics. To make its time-dependent generalization, we had to 
extend it with a real time axis. By this extension, using techniques from fibre bundles, we defined 
the biomechanical configuration bundle. On the biomechanical bundle we defined vector-fields, 
differential forms and affine connections, as well as first and second biomechanical jet manifolds 
of velocities and accelerations and prolongations of locomotion vector-fields. Using the formalism 

4 This complex geometric process is globally similar to a generic exponential decay ODE: 

Q = ~Xf(q), 

for a positive function /(<?). We can get some insight into its solution from the simple exponential decay ODE, 

q = —Xq with the solution q(t) = qoe~^, 
where q = q(t) is the observed quantity with its initial value go and A is a positive decay constant. 
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of jet manifolds of velocities and accelerations, we have developed the time-dependent Lagrangian 
biomechanics. Finally, we haver showed that the underlying geometric evolution is given by the 
parabolic Ricci flow equation. 
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